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In  the  course  of  a  lengthy  series  of  observations  since  1975.  a  large,  continuous  decrease  of  the  brightness 
of  the  solar  aureole  has  been  found  west  of  Boston,  at  Hanscom  Air  Force  Base  and  at  Lexington. 
Mass.  This  points  to  the  virtual  disappearance  from  the  lower  atmosphere  of  giant  particles  larger  than 
~  10  pm  in  size  while  total  suspended  particulates  in  Boston  and  other  U  S.  cities  have  barely 
decreased.  Results  of  calculations  to  better  understand  the  relation  between  forward  scattering  and 
aerosol  mass  distribution  [coarse  fraction  (CF)j  tire  presented.  In  addition,  a  method  to  modify  steep 
scattering  functions  calculated  for  a  plane-wave  source  iSun  treated  as  a  stan  to  those  of  the  actual  and 
limb-darkened)  Sun  is  presented.  The  calculated  wavelength  dependence  of  extinction,  which  is  lower 
than  that  observcu,  is  found  to  be  tittle  affected  by  the  CF,  but  seems,  like  forward  scattering,  to  be 
sensitive  to  mass  distribution  of  sizes  of  <  0.4  and  -06  pm  because  of  t  he  anomal  scattering  behavior  of 
spheres. 
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1.  introduction 

Over  most  of  the  globe,  the  sky  brightness  close  to  the 
Sun,  the  aureole,  is  generally  high  because  of  the 
presence  of  giant  particles  (GP),  dust  of  ~  10-200-p.m 
diameter.  Over  the  continents,  GP  from  the  ground 
or  from  industrial  activity  and  sometimes  pollen 
grains  from  forest  trees  are  carried  upwards  by 
turbulence  and  thermal  convection.  However,  the 
largest  GP,  including  pollen,  may  fall  back  to  the 
surface  during  the  night. 

Measuring  aureoles  at  angles  <  0.10°  from  the  Sun 
by  a  photometer  or  a  camera  is  difficult;  diffraction 
from  apertures  can  be  eliminated  only  by  corona- 
graphic  setups.  But  the  brightness  of  aureoles  and 
their  changes  can  be  studied  by  a  simple  method: 
Standing  at  the  very  edge  of  the  shadow  of  a  roof  at 
least  5  m  away,  one  determines  which  step  of  a  graded 
set  of  neutral-density  (ND)  filters  just  prevents  glare 
of  the  sky  at  the  Sun’s  edge.1  The  ND  value  of  this 
step  designates  the  strength  of  the  aureole.  Care 
must  be  taken  in  order  that  the  adaptation  of  the  eye 
to  moderate  outside  brightness  does  not  change. 
Related  psychological  measurements  were  discussed 
by  Holladay2  and  others.3  4  While  there  generally  is 
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no  glare  (ND  =  0)  2°-4°  from  the  Sun,  persistent 
values  of  ND  of  ~2.0  (and  sometimes  much  higher) 
have  been  observed  earlier  by  the  author  at  the  solar 
rim  at  many  locations  in  the  U.S.  and  in  Central 
Europe,  both  in  cities  and  in  the  countryside. 

2.  Decline  of  Aureole  Brightness 

The  observations  and  the  meteorological  and  environ¬ 
mental  aspects  of  the  results  obtained  during  my 
15-year  aureole  study  are  being  reported  in  some 
detail  elsewhere.5  They  will  be  discussed  here  only 
briefly  as  an  introduction  to  scattering  calculations 
that  relate  to  problems  posed  by  the  observations. 

Continuous  observations  were  begun  in  1974  at 
Phillips  Laboratory,  Hanscom  Air  Force  Base  land 
during  weekends  in  Lexington,  Mass.  I,  -  10  km  to  the 
west  of  Boston.  Excluded  from  this  study  are  much 
narrower,  generally  much  brighter,  and  fast -chang¬ 
ing  aureoles  that  are  due  to  cirrus  or  ice  crystals  from 
supercooled  water  clouds.  (Indeed,  the  method  would 
be  useful  to  detect  invisible  cirrus.)  Day-to-day  and 
diurnal  variations  of  aureole  brightness  can  be  signif¬ 
icant,  as  shown  in  Fig.  1 .  Average  annual  ND  values 
were  initially  1.5.  In  May,  pollen  caused  average 
ND’s  of  2,  and  on  windy  summer  days,  fibers  as  well 
as  stellate  hairs  from  leaves  of  forest  frees  contrib¬ 
uted  to  forward  scattering. 

The  downward  trend  of  the  ND,  in  both  back¬ 
ground  and  extreme  values,  was  realized  in  1983, 
convincing  me  to  continue  and  even  intensify  the 
observations.  Indeed,  cases  with  no  enhanced  bright- 
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Fig,  1.  Individual  aureole  observations  from  March  to  June,  for  1977  (lefts,  for  1987  ■  center,  and  for  1991  right  Ordinate  is  ND; 
brightness  goes  from  1  (no  glarei  to  1000.  Bottom  panels  are  for  morning  hours,  and  top  panels  are  for  afternoon  Heavy  dots  at  left  of 
the  panels  denote  averages  for  the  periods  November  to  January.  Data  are  for  an  aerosol  optical  thickness  (AOT  of  •-  0.25.  but  Nil's  for 
higher  turbidities  are  similar. 


ness  near  the  Sun,  i.e.,  an  apparent  lack  of  GP,  have 
become  more  and  more  frequent  in  recent  years. 

The  arrival  in  October  1991  of  the  stratospheric 
veil  of  Pinatubo  aerosol  (villi  an  AOT  of  -  0.2) 
doubled  the  brightness  of  the  sky  background  within 
10°  of  the  Sun.  Brightness  enhancement  close  to  the 
sun  was  virtually  absent,  even  during  afternoons 
(suggesting  that  convection  was  more  damped  be¬ 
cause  of  the  veil),  until  pollen  appeared  the  following 
May. 

Since  aureoles  observed  at  or  after  noontime,  when 
it  became  windy,  often  indicate  advection  from  dis¬ 
tances  of  100  to  more  than  500  km,  one  must 
conclude  that  both  day-to-day  features  and  decline 
are  not  local  phenomena  but  are  representative  of  a 
large  region. 

The  data  show  that  the  (linear)  reduction  of  aureole 
brightness  (or  of  the  vertical  load  of  GP)  between 
1977  and  1980  was  ~  35%  per  year  and  after  1980 
was  at  least  15%  annually.  In  contrast,  precipitation 
chemistry6  and  sampling  of  total  suspended  particu¬ 
lates  (TSP)  by  high-volume  filter  samplers  in  cities  in 
the  late  1980’ s  showed  reductions  in  the  eastern  part 
of  the  U.S.  of,  at  the  most,  6%  per  year7  while 
visibility  degradation8  ceased  and  TSP  in  Massachu¬ 
setts  (from  Ref.  9  and  my  evaluations  mentioned 
below)  did  not  really  change. 

However,  the  aureole  data  are  a  strong  proof  that 
GP  decreased  by  possibly  more  than  a  factor  of  100, 
which  shows,  at  least  in  part,  the  success  of  environ¬ 
mental  cleanup  efforts  in  regard  to  coarse  dust,  such 
as  that  from  power  plants  and  industry.  However, 
there  are  strong  indications  that  coarse  street  dust 
stirred  up  by  cars  on  highways  and  in  cities  also 
makes  a  significant  contribution  to  GP,  but  it  is  less 
clear  how  its  quantity  could  have  decreased. 

3.  Aerosol  Mass  Distributions  and  the  Coarse  Fraction 

The  virtual  disappearance  of  aureoles  during  my 
observations  made  it  desirable  to  try  to  relate  them  to 
other  data  that  could  provide  information  on  size 
distributions  and  possibly  on  changes  over  time. 


This  also  suggests  the  need  to  calculate  scattering 
functions. 

Indeed,  simultaneous  measurements  from  1983  to 
1986  in  several  U.S.  cities  of  TSP  and  of  inhalable 
particles  of  less  than  10  pm  in  size,  called  PM- 10, 
yielded  average  coarse  particle  fractions  ( CF’s  j,  which 
are  given  by 

CF  =  1  -  ((PM-10) /TSP),  (It 

of  0.5  to  0.6. 710  11  Such  measurements  were  also 
made,  from  1985  to  1990,  on  roofs  at  several  sites  in 
Massachusetts,  lasting  24  h  on  every  sixth  day. 
Concentrating  my  evaluation  of  printouts'*  of  the 
data  on  summer  months,  I  came  to  the  principal 
conclusions  that: 

(i)  TSP  and  the  CF  generally  were  the  same  in 
downtown  Boston  and  in  cities  located  10  to  70  km 
away.  TSP  appeared  to  be  proportional  to  turbidity 
(sun  photometry)  at  the  Phillips  Laboratory. 

(ii)  The  CF’s  of  sites  only  a  few  blocks  apart  in 
downtown  Boston  sometimes  correlated  well  (with  a 
range  of  CF’s  from  0.5  to  0.7). 

(iii)  Otherwise,  the  average  CF  was  0.6.  The  few 
data  sets  covering  the  whole  period  at  one  site  do  not 
seem  to  support  a  change  in  the  CF. 

(iv)  Farther  from  Boston,  there  are  indications  of 
a  lower  CF.  (Similar  results  were  reported  from 
Wisconsin11}.  At  the  remote  Quabbien  summit  (300 
m)  in  western  Massachusetts,  summer  values  aver¬ 
aged  to  0.4,  but  in  winter  time,  which  is  even  cleaner, 
PM-10  were  often  higher  than  TSP,  making  the  data 
suspect. 

More  detailed  mass  distributions  have  been  mea¬ 
sured  in  parallel  to  CF’s  in  some  cities;  a  typical 
average  is  the  one  shown  as  model  A  in  Fig.  2(a).10 
However,  more  sporadic  data  from  other  cities,  espe¬ 
cially  those  made  with  rotating  impactors  during  the 
day,  have  GP  peaks  shifted  to  30  or  100  p,m  (model  N 
in  Fig.  2(a) j.  I  mention  only  measurements  made  in 
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Fig.  2.  Typical  mass-size  aerosol  spectra.  Areas  under  curves 
are  proportional  to  mass.  The  (here  idealized)  10-pm-dtameter 
cutoff  defines  the  CF.  (a)  Model  A,  average  from  Phoenix  II,  is 
typical  of  U.S.  cities  listed  in  Ref.  10.  Model  N  was  obtained  by  a 
rotary  impactor  in  Seattle,  Wash.13  Model  C  was  designed  for 
large  CF.  Models  A*,  N \  and  C "  have  bigger  particles  than 
models  A,  N,  and  C,  respectively.  Model  A’  has  a  small-mode  mass 
distribution  that  produces  an  anomalous  extinction,  (b)  Models 
R,  LI,  and  L2  are  typical  for  results  at  rural  locations17  and  from 
flights.18 19 

Seattle,  Wash.,13  Chicago  and  Argonne,  Ill.,14  and  St. 
Louis,  Mo.15  Most  of  these  do  not  cover  particles 
smaller  than  4  ptm,  but  the  CF  was  most  likely  0.8  or 
larger.  The  difference  in  the  CF  can  partly  be  attrib¬ 
uted  to  sample  duration  (24  h  versus  a  few  daytime 
hours).  In  addition,  the  rain  hood  of  TSP  samplers 
causes  a  premature  deposit  of  particles  >30  pm. 
However,  all  these  data  refer  to  city  aerosols  at  busy 
street  intersections.  (In  Chicago,  more  than  of  509c 
of  the  GP  appeared  to  come  from  vehicular  street 
traffic.16)  Indeed,  for  data  from  rural  sites17  and  in 
flight,1819  the  mass  in  the  accumulation  mode  (size 
<2  pm)  and  in  the  GP  mode  was  approximately 
equal.  The  CF  is  small  because  few  GP  are  larger 
than  10  pm  [model  R  in  Fig.  2(b)]. 

For  the  translation  of  these  mass  distributions  into 
scattering  functions,  we  will,  for  simplicity  and  from  a 
lack  of  knowledge  of  vertical  distributions,  assume 
three  general  situations: 

31  Model  N  with  a  high  CF,  possibly  correspond¬ 
ing  to  strong  aureoles  measured  in  the  1970’s. 

(ii)  Model  A  with  a  CF  of  ~  0.5,  which  is  typical  of 
moderate  aureoles. 

(iii)  Flat  mass  distribution  (model  R,  rural)  with  a 
low  CF,  which  probably  corresponds  to  the  present 
aureole  data. 

Below  we  see  that  the  parameterization  of  mass 
distributions  by  CF  is  also  helpful  for  forward  scatter¬ 
ing. 


However,  it  cannot  be  overlooked  that  rural  models 
conflict  with  the  recent  CF  data  from  Massachusetts 
On  the  other  hand,  there  are  reasons,  from  only  a  few 
opportunities  of  observations,  that  aureoles  in  Boston 
also  became  weak  i  in  the  I960’s.  ND  readings  of  2  to  3 
were  common:.  We  could  possibly  solve  these  prob¬ 
lems  by  assuming  that  the  vertical  extent  of  the  cloud 
of  coarse  anthropogenic  aerosol,  over  Buston  as  well 
as  over  the  countryside,  greatly  decreased  in  recent 
years.  However,  this  would  imply  an  unlikely  clima¬ 
tological  change:  a  decrease  in  turbulence  and  con¬ 
vection  near  the  surface.  Yet  this  change  could  have 
been  subtle  because  of  the  large  fall  speed  of  the  large 
particles.  An  analogy  is  windswept  snow;  it  settles 
as  soon  as  the  gust  subsides.  A  related  experience  is 
that  aureoles  are  always  brightest  on  windy  days, 
when  soil  might  have  been  the  main  source  of  the 
dust;  but  why  then  did  they  also  fade  over  the  course 
of  the  years?  Could  it  have  to  do  with  less  open  iand 
and  less  farming? 

4.  Computation  of  Scattering  Functions  and  Results 

Calculations  of  scattering  functions,  which  use  the 
above  models  of  size  distributions,  were  made.  But 
first  a  method  is  presented  and  later  applied  to  derive 
the  effect  of  the  width  of  the  solar  disk  and  of  limb 
darkening  on  forward-scattering  calculations. 

A.  Effect  of  the  Solar  Disk  on  Scattering  Functions 

To  relate  sky  brightness  to  variations  of  the  vertical 
load  and  the  size  distribution  of  the  GP.  one  must 
consider  that  the  scattering  function  observed  near 
the  Sun  is  much  flatter  and  lower  than  if  the  Sun 
were  a  point  source.  The  closer  to  the  Sun  the 
observation  is  made  and  the  larger  the  particles  are, 
the  smaller  the  area  of  the  solar  disk  contributing  to 
the  scattering  and  the  larger  the  effect  of  limb 
darkening.  Thus  after  the  calculation  of  a  scattering 
function  P*(<t>l  (depending  on  the  scattering  angle  <i>! 
for  a  distant  point  source  (a  star:,  the  scattering 
function  must  be  converted  to  Ps< <!> 1  of  the  Sun  fthe 
angle  ik  is  measured  from  the  edge  of  the  Sun  i.  The 
function  Ps  is  expressed  as  the  convolution  of  P*  and 
the  angular  brightness  distribution  of  the  Sun. 

An  analytical  solution  of  the  integral  exists-"  for 
scattering  functions  of  the  form 

P* (<t>)  -  c  (<f>-  +  d i;  (21 

results  are  given  in  Ref.  2.  For  arbitrary  scattering 
functions,  linearization  has  been  proposed.23  Sample 
results  of  strict  solutions,  including  limb  darkening 
for  relatively  flat  forward  scattering,  are  presented  in 
Refs.  23  and  24.  However,  a  field  of  view  of  even  less 
than  0.5‘J  of  a  photometer  may  dominate  the  relation 
between  P*  and  Ps  to  up  to  a  few  degrees  from  the 
Sun.24  This  problem  is  irrelevant  to  our  method  of 
observation. 

In  this  paper,  simple  algorithms  derived  from  a 
graphic  presentation  of  the  scattering  geometry  and 
of  the  brightness  distribution  of  the  solar  disk  are 
used  for  the  convolution.  The  disk  is  divided  into 
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Fig.  3.  Subdivision  of  solar  disk  only  partly  shown  ■  to  modify  the 
small-angle  scattering  function,  which  is  caused  by  a  distant  point 
source,  to  that  from  the  limb-darkened  solar  disk.  :ai  For  ob¬ 
served  point  at  edge  of  solar  disk  'ill  =  O  i,  with  10  sectors;  b  for 
the  sky  0.2s  to  1.6°  from  Sun.  Table  1  displays  the  weights  for  the 
sectors  and  the  computation  scheme. 


ring  zones  that  are  concentric  to  the  respective 
observation  points  in  the  sky.  The  first  evaluation, 
which  is  depicted  in  Fig.  3(a),  is  for  the  sky  right  at  the 
edge  of  the  Sun,  where  the  highest  accuracy  is  needed 
if  the  scattering  function  is  steep.  Considering  that 
P*  is  calculated  in  0.1°  increments  of  <{>  (including  a 
value  at  0.04°),  P*  for  an  effective  angle  of  6  =  0.02°  is 
obtained  from  |P*( 0°)  P*(0.04°)j1  2,  and  P*  for  0.07° 
and  0.13'  are  similarly  obtained  from  the  next  pairs  of 
primary  P  intensities.  After  taking  limb  darkening 
into  consideration,  we  graphically  determine  and 
normalize  (as  for  the  coefficients  CD  in  Table  1)  the 
areas  of  the  zones  by  the  respective  sums.  The 
scheme  for  ib  =  0.04°  uses  the  same  coefficients  CD  as 
for  ib  =  0°.  Simpler,  but  basically  similar,  are  the 
evaluations  for  tb  =  0.1°,  with  coefficients  CE.  The 


scheme  for  <b  =  0  2',  which  uses  angles  <b  to  0.7  and 
yields  coefficients  CC,  is  shown  in  Fig.  3sb  The 
same  scheme  is  repeated  for  starting  angles  of  ib  = 
0.3\  0.4  ■  •  •  1.6',  so  that  angular  increments  of  6  can 
be  sand  are  here;  larger  beyond  2. 1  . 

Although  the  errors  of  the  coefficients  Table  1 
might  be  as  much  as  30,  this  is  expected  to  cancel  in 
the  calculations.  Indeed,  for  a  P*  of  constant  value, 
P  is  fiat  to  within  0.50.  Results  of  the  method  for 
model  A  can  be  seen  in  Fig.  4.  The  flattening,  which 
is  -400  at  the  Sun’s  edge,  obviously  would  be 
detrimental  to  attempts  to  recover  by  inversion  the 
true  scattering  function  For  particles  such  as  ice 
crystals  with  a  relatively  narrow  size  distribution 
peaking  at  a  diameter  D  of  220  pm.  the  edge 
brightness  is  reduced  much  more,  to  20  only  !  not 
shownc 

The  validity  of  the  scheme  extends  to  size  parame¬ 
ters  2t ~f  *  D  X  of  at  least  1800.  This  follows  from 
Eq.  i  3a i  of  Table  1  in  that  the  intensity  in  the  central 
diffraction  peak  should  differ  little  at  d>  =  0.01s  from 
that  at  <b  =  0.0°. 

The  discontinuity  seen  in  Fig.  4  at  6  =  1.6'  has  a 
simple  explanation:  the  last  modified  intensity  is 
derived  from  an  effective  angle  of  ib  =  1.8  .  and 
therefore  is  lower  than  the  next  unmodified  value  at 
ib  =  1.7°.  The  jump  is  largest  where  the  primary' 
scattering  function  is  steepest !  solid  curve  ).  Carrying 
the  smoothing  on  to  larger  angles  will  eliminate  the 
jump. 

B.  Some  Scattering  Models 

Calculations  of  scattering  functions  that  use  the 
above  models  of  size  distributions  were  made  for 
particle  sizes  from  0.16  to  540  pm  in  diameter  by  a 
reliable  program,  although  calculations  for  size  param- 


Table  1 .  Setup  of  Calculations  to  Convert  Scattering  Functions  for  a  Starlike  Sun  ( P*  i  to  those  of  the  Acfual  Sun  t P$  < 


Primary  angles  6  for  which  P* l/v ),  henceforth  called  PuV !,  is  assumed  to  have  been  calculated 
0.0,0.04,  0.1,  (0.1  >  •  to  2.1  deg,  (Pill  to Pi21|j 


Auxiliary  scattering  intensities  and  approximate  angles: 


P02  =  [miPi2)]>  2  4  =  0.02°. 

i3a> 

P07  =  ]Pi2lPl3)|!  2  4  =  0.07°, 

(3bi 

P13  =  [P(3i  P(  4 )  ] 1  2  4  =  0.13°, 

i-3c 

1*11)  =  CD1  P(l)  +  CD2P02  +  CD3  P(2i  +  CD4  P07  +  CD5P(3i 

+  CD6  P13  +  CD7  P(4l  +  CD8  Pi 5)  +  CD9  P(6>  +  CD10  PI 7i, 

O 

o 

o 

II 

-3 

i  4a  i 

Ps( 2 1  =  CD1  P(2 1  +  CD2P07  +  CD3  P07  +  CD4P13)  +  CD5  P13 
+  CD6  P( 4 i  +  (CD7  +  CD8)  Pi 5)  +  CD9PI6)  +  CD10PI7). 

4  =  0.04" 

i4b( 

PSI3)  =  CE1  P(3)  +  CE2  P13  +  CE3  P(4)  +  CE4  SQRT(P<4)  P(5tj 
+  CE5  P(5|  +  CE6  P<6)  +  CE7  PI 7)  +  CE8  P(8l. 

4  =  0.10° 

1 4c  i 

P'i(/)  for  /=  4  to  /  =  18: 

Ps!/l  =  CC1  Pi/)  +  CC2  Pil  +  1)  +  CC3  Pi  /  +  2) 

+  CC4  P(/  +  3)  +  CC5  Pi/  +  4)  +  CC6  P(/  +  5) 

4  =  0.2°  to  1.6°. 

{ 5 1 

Coefficients  CD  |for  Ps  at  ill  =  0.0°  and  0.04°,  Eq.  i4a)|: 

0.0013.  0.072,  0  120,  0.023, 0.C52,  C.119,  0.102.  0.312,  O.ktw.  0.043 
Coefficients  CF,  jforPs  at  >1<  =  0.1,  Eq.  (4b)]: 

0.0083,  0.033,  0.076,  0.130,  0.210,  0.290,  0.217,  0.030 
Coefficients  CC  |  for  Ps  at  ill  =  0.2  to  1 .6,  Eq.  (4c)]: 

0.052,  0.192,  0.281,  0.281,  0.252,  0  032 
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Fig.  4.  Comprehensive  presentation  of  scattering  data.  Shown 
are  scattering  functions  for  model  A  at  X  0.5  pm  (P*,  unmodified, 
and  Ps  modified  by  the  solar-disk  scheme),  and  P5  for  X  0.86  and  1 .6 
urn.  The  bottom  part  of  the  graph  shows  color  ratios.  The 
angles  had  to  be  increased  by  0.02°  for  plotting.  EPA  is  the 
effective  extinction  cross  section.  X  exponents  for  the  X  pairs  0.86, 
0.50  and  1.6,  0.86  were  0.69  and  0.58,  respectively. 

eters  >  500  and  angles  <  5°  were  usually  replaced  by 
the  corresponding  Bessel  functions  from  Fraunhofer 
diffraction  optics.  The  logarithmic  increment  (base 
10)  in  size  is  0.1.  Initially,  a  power-law  size  distribu¬ 
tion  is  assumed  (the  label  EXP  in  Fig.  4),  but  specify¬ 
ing  factors  in  logarithmic  size  increments  of  0.26  will 
result  in  a  logarithmically  interpolated  mass  distribu¬ 
tion.  Calculations  were  usually  restricted  to  the 
three  wavelengths  (of  Fig.  4,  but  in  the  following, 
generally  only  results  for  K  =  0.5  pm  are  considered. 
As  is  appropriate  for  the  sky  near  the  Sun  in  typical 
situations  of  continental  pollution,  no  molecular  scat¬ 
tering  was  included.  By  remaining  in  the  realm  of 
normalized  scattering  functions,  actual  mass  concen¬ 
trations  do  not  matter,  and  effects  that  are  due  to 
variations  of  the  path  length  in  the  atmosphere  (solar 
elevation  angle)  and  aerosol  extinction  are  neglected, 
as  is  multiple  scattering.  In  short,  the  results  are 
generally  valid  for  small  overall  turbidity  and  the  Sun 
near  zenith. 

The  results  of  the  calculations  for  the  edge  of  the 
Sun,  plotted  against  the  CF,  are  shown  in  Fig.  5. 
On  the  right-hand  side,  an  ordinate  relating  P  to  sky 
radiance  B  with  respect  to  extraterrestrial  radiance 
S0  of  the  solar  disk  (of  solid  angle  of  67  x  10~6  sr),  has 
been  added  (e.g.,  Refs.  23  and  24).  A  typical  AOT  of 
0.1  at  X  =  0.5  pm  and  an  air  mass  of  1.5  are  assumed. 
The  ND  scale  is  based  on  measurements  taken  in  a 


Fig.  5.  Forward  scattering  versus  CF.  Normalized  forward  scat¬ 
tering  P*  and  P1'  at  the  edge  of  the  Sun  for  some  mass  distribution 
models  of  Fig.  2  are  shown.  The  right-hand  side  shows  ordinates 
for  the  radiance  ratio  of  sky  to  Sun  and  for  the  ND.  w  ith  an  AOT  of 
0.1  and  an  air  massofl.5.  The  upper  ends  of  the  data  bars  are  for 
a  starlike  sun.  and  the  lower  ends  are  for  the  actual  Sun  Dashed 
curves  indicate  the  downshift  of  several  models  caused  by  mass 
distribution  A'  of  their  small  mode.  At  CF  ■  0.2.  P  drops  off 
steeply  to  low  scattering  by  power-law  size  distributions  and  to 
values  marked  as  G  and  V.  which  was  observed  in  the  Alps-'  ,r.  the 
expected  Rayleigh  scattering  t  RAY i  is  indicated. 


clear  sky  with  a  high  Sun:  ND  =  0  corresponds  to 
B/S0  -  40  xg  10~6  ±  20#. 

Data  points  obtained  from  models  A.  N,  and  C.  and 
with  an  increase  in  the  height  or  width  of  the  GP 
mode  (models  A",  N+,  and  C ’ ),  together  with  models 
LI,  L2,  and  R  from  Fig.  2(b),  are  plotted  against  CF. 
P*  and  Ps  are  shown  as  upper  and  lower  ends  of  the 
data  bars,  respectively.  For  CF  >  0.3,  brightness  is 
-40%  lower  for  the  actual  Sun  than  for  a  starlike 
sun. 

At  a  small  CF,  one  must  consider  that  scattering  at 
0°  of  power-law  size  distributions  with  the  exponent 
-3  (which  would  produce  a  straight  line  in  Fig,  2)  or 
flatter  distributions  would  rise  unlimited  unless  a 
largest  size  is  specified.  For  the  power-law  size 
distributions  indicated  in  Fig.  5  (mass  decreasing 
with  increasing  size),  the  CF  is  essentially  zero. 

The  result  is  that  log  P  or  NDs  an  increase,  which  is 
proportional  to  the  CF.  Between  CF  -  0.2  and  0.9, 
the  slope  d  ND/rf  CF  is  -  2.2. 

It  was  also  investigated  how  scattering  by  a  sharp 
pollen  peak  fares  if  the  GP  part  is  subjected  to  a 
uniform  reduction  (from  mass  distribution  Ql  to  Q2 
in  Fig.  6),  say,  by  fallout.  This  caused  the  CF  to 
halve,  and  decreased  the  brightness  at  angles  <  4°,  as 
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Fig.  6.  Sensitivity  study:  Effect  on  scattering  function  Ps  i\  = 
0  50  |im,  inset*  of  the  reduction  of  a  sharp  GP  peak  i  model  Ql )  by  a 
factor  of  0.25  (model  Q2l.  The  brigh'ness  at  the  Sun's  edge  is 
affected  by  the  same  factor.  The  color  ,atio  (A,  0.86  0.50lofQ2is 
slightly  lower  than  for  Ql. 


Fig.  i.  Wavelength  exponents  a  of  extinction  from  regular  models 
(in  the  figure  in  the  region  below  u  *  1.0)  and  of  A  -type  models  in 
the  upper  part  The  solid  circles  indicate  a  between  wavelengths 
of  0.50  and  0.86  am  and  the  long  dashes  indicate  a  between 
wavelengths  of  0.86  and  ’  5  um  Observations  of  the  exponent  a 
deviate  little  from  1.3  and  were  found  to  be  insensitive  to  aureole 
brightness.  Models  S  (not  mentioned  abovei  and  N  do  not  have 
the  usual  minimum  mass  near  sizes  of  -  1  nm;  the  model  02  has, 
compared  with  model  01,  a  bulge  at  D  =  0.6  p.m. 


shown  in  the  inset  of  Fig.  6.  But  the  reduction  of  the 
brightness  reflects  the  change  of  mass  at  the  peak 
which  attests  to  the  dominance  of  scattering  by  the 
pollen.  However,  the  color  ratio  of  model  Q2,  which 
also  peaks  at  1°,  is  slightly  lower  than  that  of  Ql. 
Otherwise,  the  sharp  peak  reduces  the  absolute  scat¬ 
tering  considerably  with  regard  to  CF  (Fig.  5). 

With  model  A,  which  was  modified  by  a  high 
number  of  particles  at  D  <  0.4  pm  ( model  A' ),  a  large 
increase  of  the  X  dependency  of  extinction  and  a  loss 
of  scattering  at  angles  <  30°  was  found.  Therefore, 
more  models  (C  and  C*)  were  modified,  too,  at  X  <  0.4 
pm  by  the  same  distribution  (models  C'  and  C*’). 
This  reduces  scattering  even  at  a  high  CF,  as  indi¬ 
cated  by  the  dashed  curves  in  Fig.  5.  The  effects  are 
related  to  the  anomalous  extinction  by  these  small 
particles. 

Overall,  these  results  show  that  neglect  in  our 
database  of  the  shape  and  width  of  observed  aureoles 
is  not  detrimental  if  particle  mass  is  of  prime  interest. 

Now,  what  is  the  relation  between  the  model 
results  at  the  Sun’s  edge  and  the  ND  scale  on  the 
right-hand  side  of  Fig.  5?  The  agreement  appears  to 
'  be  quite  good.  These  attempts  to  get  confirmation  of 
the  aureole  data  from  scattering  theory  seem  to 
indicate  that  the  CF  concept  is,  in  general,  quite 
helpful  for  the  study  of  forward  scattering  by  typical 
pollution  aerosols. 

5.  Wavelength  Dependence  of  Extinction 

A  brief  look  at  the  possibility  of  connecting  routine 
Sun  photometer  results  concerning  the  wavelength 
exponent  a  of  aerosol  extinction  to  aureoles  is  now 
given.  This  value  (as  usual,  its  absolute  value)  has 
the  advantage  in  sampling  the  vertical  load,  as  does 
the  ND,  and  in  referring  to  the  place  of  the  ND 
observation.  The  respective  results  for  a  of  our 
regular  models  decrease,  as  seen  in  the  lower  part  of 
Fig.‘7,  between  wavelengths  of  0.50  and  1.6  (im  from 
about  0.9  to  0.7,  as  the  CF  goes  from  0.2  to  0.9.  The 
.  -decrease  of  a  to  0  (as  during  dust  storms  in  the  visible 
part  of  the  spectrum)  occurs  only  at  much  longer 
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wavelengths.  However,  anomal  extinction  that  is 
due  to  a  large  number  of  particles  with  D  <  0.4  pm 
(as  in  model  A’  1  increases  a  for  the  considered  CF’s  to 
values  greater  than  one,  fitting  them  into  the  upper 
part  of  Fig.  7. 

The  actual  influence  of  GP  or  small  particles  on 
extinction  was  studied  by  sun  photometer  data  at 
wavelengths  of  0.38,  0.50,  and  0.86  pm  for  the 
summers  from  1976  to  1978,  in  the  absence  of 
volcanic  turbidity.  The  average  a  being  1.3,  even 
observations  with  ND  >  2.7  show  no  smaller  expo¬ 
nents  than  those  with  ND  <  1.0,  regardless  of 
turbidity.  Since  the  mass  distributions  causing 
anomal  extinction  cannot  always  prevail,  this  method 
also  does  not  help  to  establish  trends  of  ND.  Is  it 
that  the  mass  distributions  assumed  are  too  struc¬ 
tured  to  affect  a?  The  problems  also  probably  can¬ 
not  be  resolved  by  assuming  that  the  aureole  decline 
is  reflecting  the  decrease  of  particles  with  D  >  30  pm, 
which  are  known  to  have  been  missed  in  measure¬ 
ments  of  TSP  during  windy  weather. 

6.  Summary  and  Conclusions 

For  the  past  sixteen  years  I  have  estimated,  both  at 
Hanscom  Air  Force  Base  and  in  Lexington,  to  the 
west  of  Boston,  the  bright  ness  of  the  solar  aureole  at 
the  edge  of  the  Sun.  My  purpose  was  to  study 
variations  of  the  vertical  load  of  giant  aerosol  particles. 
Unexpectedly,  the  average  brightness  of  aureoles 
decreased  continually  until  they  essentially  disap¬ 
peared  in  the  background  of  small  particle  scattering. 
However,  the  interrelation  between  the  decline  of  the 
aureole  brightness,  for  which  there  is  no  support 
from  other  evidence,  and  comparatively  constant 
roof-level  pollution  is  still  not  clear.  To  partly  ad¬ 
dress  this  question,  emphasis  in  this  study  was  on 
calculations  of  forward  scattering  to  better  under¬ 
stand  the  relation  of  aureole  brightness  to  typical 
mass-size  distributions  of  aerosols.  The  calcula¬ 
tions  include  the  presentation  and  use  of  a  method  to 
account  for  the  effect  of  the  width  of  the  solar  disk 
and  of  limb  darkening  on  forward  scattering.  These 


effects  result  in  flattening  and  lowering  of  scattering 
near  the  Sun,  but  they  weaken  the  relationship 
between  the  shape  of  the  scattering  function  and  the 
size  distribution  of  particles.  Thus  it  makes  observa¬ 
tions  easier  but  would  greatly  hinder  efforts  to  invert 
measured  scattering  functions  to  yield  size-mass 
distributions. 

Several  aerosol  models  have  been  investigated. 
Surprisingly,  the  concept  of  the  CF  (mass  of  particles 
>  10  (xm/total  mass)  often  used  (in  the  form  1  -  CF) 
to  characterize  the  results  of  routine  monitoring  by 
filter  samplers  allows  one  to  present  normalized 
scattering  functions  at  the  solar  rim  in  a  seemingly 
concise  manner.  That  is,  the  CF  and  the  glare  (ND) 
of  the  aureole  are  closely  related  and  are  rather 
independent  of  the  shape  and  the  width  of  the  mass 
distribution  of  the  large  particles.  Hence  the  angu¬ 
lar  width,  shape,  and  coloration  of  the  aureole  (which 
were  neglected  in  the  presentation  of  the  observa¬ 
tions)  also  are  not  important. 

The  models  clearly  show  that  the  recent  aureoles 
observed  at  Phillips  Laboratory  (and  probably  in 
Boston  as  well)  are  caused  by  low-CF  mass  distribu¬ 
tions.  But  in  1990  roof-level  values  of  CF  ~  0.5 
were  still  measured  in  Boston,  and,  more  relevant  to 
aureoles  at  Phillips  Laboratory,  in  cities  of  central 
and  western  Massachusetts.  A  way  out  of  this  di¬ 
lemma  is  the  assumption  that  the  GP’s  are  now 
falling  out  faster  because  of  lower  turbulence  near  the 
surface.  And  should  Boston  aureoles  now  indeed  be 
as  low  as  those  at  Phillips  Laboratory,  one  might  have 
to  assume  that  the  street-corner  pollution  is  not 
representative  of  the  optics  of  the  Boston  sky,  or  even 
of  the  whole  of  the  Boston  roof-level  aerosol. 

An  excursion  into  results  on  the  wavelength  depen¬ 
dence  of  extinction  reveals  that  the  X  exponent  of  the 
regular  models  is  only  slightly  affected  by  the  CF  (and 
the  ND),  but  is  much  lower  than  observed.  Hence  it 
is  obvious  that  several  questions  remain  regarding 
the  interplay  between  mass  distributions,  optical 
observations,  and  modeling. 
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